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ABSTRACT 

Context. Radio observations are an effective tool to discover particle acceleration regions in colliding-wind binaries, through detec¬ 
tion of synchrotron radiation; these regions are natural laboratories for the study of relativistic particles. Wind-collision region (WCR) 
models can reproduce the radio continuum spectra of massive binaries that contain both thermal and non-thermal radio emission; 
however, key constraints for models come from high-resolution imaging. Only five WCRs have been resolved to date at radio fre¬ 
quencies at milliarcsec (mas) angular scales. The source HD 93129A, prototype of the very few known 021 stars, is a promising target 
for study; recently, a second massive, early-type star about 50 mas away was discovered, and a non-thermal radio source detected in 
the region. Preliminary long-baseline array data suggest that a significant fraction of the radio emission from the system comes from 
a putative WCR. 

Aims. We sought evidence that HD 93129Aisa massive binary system with colliding stellar winds that produce non-thermal radiation, 
through spatially resolved images of the radio emitting regions. 

Methods. We completed observations with the Australian Long Baseline Array (LBA) to resolve the system at mas angular resolutions 
and reduced archival Australia Telescope Compact Array (ATCA) data to derive the total radio emission. We also compiled optical 
astrometric data of the system in a homogeneous way. We reduced historical Hubble Space Telescope data and obtained absolute and 
relative astrometry with milliarcsec accuracy. 

Results. 

Conclusions. The astrometric analysis leads us to conclude that the two stars in HD 93129A form a gravitationally bound system. The 
LBA data reveal an extended arc-shaped non-thermal source between the two stars, indicative of a WCR. The wind momentum-rate 
ratio of the two stellar winds is estimated. The ATCA data show a point source with a change in flux level between 2003-4 and 
2008-9, that is modeled with a non-thermal power-law spectrum with spectral indices of -1.03 ± 0.09 and -1.21 ± 0.03 respectively. 
The mass-loss rates derived from the deduced thermal radio emission and from the characteristics of the WCR are consistent with 
estimates derived by other authors. 

Key words. Stars: massive - binaries: general - Radio continuum: stars - Radiation mechanisms: non-thermal - stars: individual: 
HD 93129A 


1. Introduction 

Wind-collision regions (WCR) in binary systems with two hot 
massive stars constitute natural laboratories to study relativis¬ 
tic particles ( Eichler & Usov||1993t [Benaglia & Romero|[2003 


|Pittard & Dougherty||2006| |De Becker||2007| l, as do supernova 

remnants but with somewhat different physical and geometri¬ 
cal parameters. A WCR is the site of particle acceleration that 
generates a population of relativistic electrons, as now demon¬ 
strated in more than 40 massive binary systems. The particle- 
accelerating colliding-wind binary (PACWB) status of these ob¬ 
jects is in most cases revealed by radio observations through the 
detection of synchrotron emission (|De Becker & Raucq|2013|l. 


WCR models are able to reproduce the radio continuum 
spectra of these massive binaries ( |Pittard & Dougherty||2006] l. 
The radiation from these systems is the superposition of ther¬ 
mal emission from the stellar winds of each binary component 
and the non-thermal emission from the WCR. A key ingredi¬ 
ent for constraining models of the non-thermal component is ex¬ 
pected to come from high-resolution radio imaging. The non- 
thermal emission region associated with colliding winds has 
been mapped at radio frequencies for WR 140 (|Dougherty et 
ar||2005|), W R 146 ( [O’Connor et~ari|2005|l WR 147 ([Williams 


et al.[[1997[ l, CygOB2#5 ( [Contreras et aTl[T997 


al.|2()l l[ l~and Cyg OB2 #9 ( [Dougherty & Pittard 2006| l. In each 


Ortiz-Leon et 


case, the WCR is resolved by observations with angular resolu¬ 
tions of a few milliarcseconds (mas). 
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HD 93129A is the brightest source in the most crowded part 
of Tr 14 in the Carina nebula and has been studied for several 
decades. The record of its special spectral characteristics can be 
traced back to [Payne | ( |1927| ) (see also Table [T] for information 
on other studies). It now is recognised as hosting the prototype 
02 If+ ( jWalborn et al.|2002l l. Radio emission from HD 93129A 
has been detected between 1.4 and 24.5 GHz (e.g. [Benaglia et| 
|ar]|2006l ), showing a flux decreasing with frequency indicative 
of synchrotron radiation, most probably produced in a WCR be¬ 
tween stars in a colliding-wind binary system. The Australian 
Long Baseline Array (LBA^observation presented by [Benaglia 
et al.| pOlO I suggested that a significant fraction of the radio 


emission comes from a region compatible with a WCR. 

In this paper, we offer an in-depth discussion of a more recent 
LBA observation of HD 93129A, aiming to provide compelling 
evidence of the WCR origin for the bulk of the radio emission, 
in agreement with the so-called standard scenario for PACWBs. 
These results allow us to spatially resolve for the first time the 
synchrotron emission region and to derive specific information 
about the non-thermal radio contribution. The two component 
nature of HD 93129A is presented, including the most recent as¬ 
trometric measurements that show the relative position of the 
components over time (Sect.j^, and a new analysis of Hubble 
Space Telescope/Fine Guidance Sensor (HST/FGS) data that al¬ 
low us to correlate the optical positions of the two components 
with the position of the radio emission (Sect.j^. The LBA ob¬ 
servations are presented in Sect.® The main results and general 
discussion are in Sections |^and^ with a final statement of our 
conclusions in Sect. [7] 


2. The HD93129A stellar system 


The object HD93129A was first classified as an 03 star (Wal- 
[bom|1987]l, though later was reclassified as the prototype 02 I 
star |Walborn et ^ 2002| l. Mason et al. ( 1998[ l had suggested 
it was a speckle binary, and was confirmed as having two com¬ 
ponents, Aa and Ab, from HST/FGS observations ( [Nelan et al.| 
2004 [ l; these authors provide the separation at the epoch of the 
observations, the difference in apparent visual magnitude be¬ 
tween the primary and the secondary and the position angle (PA) 
(see also Nelan et al. [2010[ l. Maiz-Apellaniz et al. ( 2008| l de¬ 
tected proper motion along the radius vector between the compo¬ 
nents after analysis of data from 1996, 2002 and 2004; the binary 
nature is supported by the relative velocity of about 2 mas yr“*, 
higher than the expected relative velocities in the held. The au¬ 
thors suggested a highly elliptical or inclined orbit, or both, al¬ 
though they did not determine the orbital parameters. The com¬ 
panion is assumed coeval with the primary and of the same spec¬ 
tral type as various other OB stars nearby (03.5 V) (s ee TablefTjl. 
Through detailed ultra-violet and Ha line modeling, Taresch et 
(aL[([1997[| and[Repolust et al.[([2004[) derived the main parameters 
of HD 93129A (see Table [IJ). 

The relative motion of the two stars can be traced by archived 
HST/FGS and HST/ACS0 observations from 1996 to 2009, 
along with recent VLT/NACO and VLTI/PI0NIER[^ observa¬ 
tions of Sana et al. (2014). Figure[^shows the relative astrome¬ 
try from these measurements over time. The position angles for 


* The Australian Long Baseline Array is part of the Australia Tele¬ 
scope which is funded by the Commonwealth of Australia for operation 
as a National Facility managed by CSIRO. 

^ Hubble Space Telescope/Advanced Camera for Surveys. 

^ Very Large Telescope NAos-COnica and Precision Integrated-Optics 
Near-infrared Imaging ExpeRiment. 
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1996 and 2002 have large uncertainties because the binary com¬ 
ponents were resolved along only one FGS axis, due to the HS T 
roll angle resulting in a small projected separation along the un¬ 
resolved axis. For later epochs, the approximate position angle 
was known, so the observations were planned to avoid this prob¬ 
lem. 

HD 93129A was included in the recent Galactic O star Spec¬ 
troscopic Survey of Sota et al. ([2014| l. The coordinates from the 
Tycho catalogue ( Hog et al.|2()00| l assume a single star since the 
presence of a companion star was unknown at the time. Hippar- 
cos did not resolved the system. These coordinates have errors 
of about 0.1”, an order of magnitude larger than the expected 
angular resolution of the LBA data of ~ 10 mas. 


Table 1. Adopted parameters of HD93129A, relevant here. 


Parameter 

Value 

Unit 

Ref. 

Aa Spectral type 

02 If-H 


Wal02 

Ab Spectral type 

03.5 V 


Wal06 

Am 

0.9+0.05 

mag 

NellO 

Wind terminal velocity t 

3200+200 

km s“' 

Tar97 

System mass 

200 + 45 

Mo 

Mai08 

Effective temperaturef 

42500 

K 

Rep04 

System luminosity (log) 

6.2 

L/Lq 

Rep04 

X-ray luminosity ft 

1.3 X 10-^ 

^bol 

Cohll 

Distance 

2.5 

kpc 

Wal95 


f: Derived if a single star, ft: Derived from the energy 
band (0.5 - 8) keV. 


References. Wal02: 

Walbom et al. ( 2002); Wal06: Walborn, priv. com- 

munication; NellO: 

Nelan et al.H2010|l; Tar97: [Taresch et al. 

19971; 


Mai08: [Maiz-Apellaniz et al.[ ( [2608lirRep04: [Repolust et al. 


Coh 11: [Cohen et al.[ ( |2011| >; Wal95: |Walbom| ( fl995r 


120041; 


3. Hubble Space Telescope/FGS data analysis 

3.1. Astrometry 

To improve upon the coordinates of HD 93129A from the Tycho 
catalogue, we re-analyzed HST/FGS measurements. The three 
Fine Guidance Sensors on the HS T are 2-channel (x,y) white- 
light shearing interferometers that provide line-of-sight posi¬ 
tional measurements of guide stars to enable accurate pointing 
and stabilization of the HST. The FGS is capable of resolv¬ 
ing sources down to ~15 mas when operated in its high angular 
TRANS (or Transfer) mode and sub-millarcsecond astrometry 
when operated in POS (or Position) mode. A description of the 
FGS and its operation can be found in [Nelan[ ( [2014 l. 

In eight epochs between 2006 and 2009, FGS observed 
HD93129A in TRANS mode, along with “wide angle” POS 
mode observations of the binary itself and seven reference held 
stars within one arcminute (see Fig. [^. The TRANS mode ob¬ 
servations resolve HD 93129Aab, providing the component sep¬ 
aration, position angle, and magnitude difference at each epoch. 
The POS mode observations supply the relative positions of the 
held stars and composite HD 93129A system. Combining data 
from these two observing modes gives the position of Aa and Ab 
relative to these held stars at each epoch. Combining all epochs 
yields a local FGS catalog of the relative positions and proper 
motions of reference stars and binary components. 

Figure shows the interference fringes of HD 93129A for 
the Jan 27, 2009 observations, the best fitting binary star model 
and, for comparison, the fringes of an unresolved point source. 
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Table 2. PPMXL and FGS astrometry of the HD 93129A field. 


Star 


RA(h) 


PPMXL 
Dec (deg) 




Us 


FGS 

T] 


l^a 


IJ^S 


R1 10.7325347 (13) 

R2* 10.7322835 (14) 

R3 10.7340882 (13) 

R4 10.7344007 (49) 

R6 10.7326563 (13) 

R7 10.7315615 (13) 

R8 10.7296359 (15) 


-59.535576 (19) 
-59.541006 (16) 
-59.541208 (19) 
-59.552781 (74) 
-59.560735 (19) 
-59.557925 (19) 
-59.548569 (15) 


-6.6 (2.5) -0.6 (2.5) 


-4.8 (3.2) 
24.6 (17) 
-4.8 (2.1) 
-6.6 (2.5) 
-6.4 (2.0) 


0.6 (3.2) 

60.5 (17) 
3.0(2.!) 
1.6 (2.5) 

2.5 (2.0) 


-2.6531 (0.9) 
-9.5351 (1.0) 
39.9186 (1.2) 
48.0711 (1.1) 
0.8902(1.2) 
-29.0924(1.1) 
-81.8936 (1.6) 


43.2955 (0.5) 
23.4780 (0.5) 
23.1199 (0.5) 
-19.2493 (0.5) 
-47.2647 (0.4) 
-37.1969 (0.2) 
-3.6759 (0.3) 


-5.78 (0.29) 
-5.78 (0.53) 
-6.73 (0.32) 
-5.34 (0.44) 
-6.29 (0.32) 
-5.46 (0.29) 
-5.53 (0.32) 


1.54 (0.31) 
1.71 (0.50) 
1.38 (0.36) 
2.22 (0.45) 
1.64(0.34) 
1.60 (0.30) 
2.12(0.30) 


‘Star’; Reference star identification, same as in Figure|^ All errors are shown within parenthesis. PPMXL coordinates 

(J2000): RA in hours, Dec in degrees, with errors in mas. The J2000 FGS tangential plane coordinates, rj, are given in 
arcsec, relative to the point mid-way between the binary components; errors are in mas. Proper motions (ju) in mas yr ’, 
with errors in the same unit. *; Reference star 2 is not in the PPMXL catalog, so coordinates from UCAC4 are quoted. 
The FGS proper motions are based upon FGS data but are constrained by the PPMXL values. The entries in this table 
serve as the input values for the tangent plane analysis outlined in Section]^ 



Fig. 1. The relative positions, including error bars, of the two compo¬ 
nents of HD93129A, re-derived from archive HST/FGS (blue labels) 
data, along with VLT/NACO and VLTI/PIONIER (red labels: taken 
from 2011.18 to 2013.09) data ( [Sana et ^[2014^ . The epoch of the 
measurement is given for each point. 


The binary model takes into account the light contributed by 
HD93129B which, 2.73” away, accounts for about 5 percent of 
the light of HD93129A (FGS optics collimates a source onto a 
circular beam with a full width-half maximum intensity diameter 
of about 3”). Figurej^also shows the location on the fringes that 
the POS mode observations of HD93129A “lock” on to. The 
model fit to the TRANS mode observations provides the offset 
of each component from this point, which thereby allows us to 
create pseudo POS mode exposures for both HD 93129Aa and 
HD 93129Ah, and hence accurate positions for each star relative 
to the reference stars. The bright O star HD 202124 was used as 
the point-source calibrator since it is nearly the same B-V color 
as HD 93129A, and it was observed by FGS only eight days after 
HD 93129A. The orientation of HST at the time of the observa¬ 



Fig. 2. The core of Tr 14 and the selected FGS reference stars used in 
the astrometric measurements. The image (in log stretch) is from HST 
proposal 10602 using ACS/WFC with F550M, which is similar to the 
central pass band of the FGS measurements. The two stars HD93129A 
and HD93129B, separated by 2.7”, are not clearly resolved in this im¬ 
age since the roll angle of HS T at the time caused their position an¬ 
gle to be aligned with the diffraction spike. In this image the compo¬ 
nents of HD 93129A are not resolved. Note that the separation between 
HD93129A and HD 93128 is about 24”. North is up and East is to the 
left. 


tion caused the projected separation of HD 93129A and B along 
the FGS y-axis to be only 0.233", as is evident in the figure. 

To compare the FGS coordinates of HD 93129Aab to the ra¬ 
dio source, it is necessary to convert the relative FGS positions 
to absolute International Celestial Reference System (ICRS) co¬ 
ordinates. The FGS coordinate frame lies in a tangential plane, 
with the zero point arbitrarily set to the point mid way between 
the binary components (at the zero point the FGS Cartesian plane 
is ta ngential to the celes tial sphere). We used the PPMXL cata¬ 
log ( RDeseTetarj^O^ for the coordinates and proper motions 
of the reference stars as input to a model that finds the RA,Dec of 
this tangent point. This is a two step process. The FGS observa¬ 
tions were obtained between 2006 and 2009, while the PPXML 


PPMXL is a catalog of positions, proper motions, 2MASS and opti¬ 
cal photometry of 900 million stars and galaxies. 
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Fig. 3. Interference fringes of HD 93129A for Jan 27, 2009 FGS obser¬ 
vations (solid line), the best fitting binary star model (dashed line) and, 
for comparison, the fringes of an unresolved point source (dotted line). 
The "X" shows the location of the composite fringe that POS mode ob¬ 
servations tracked. Note that due to the roll of HST at the time of the 
observations, HD93129B has a projected separation from A along the 
y-axis of only 0.233", although its true distance is 2.73". 


coordinates are epoch J2000. The proper motions of the refer¬ 
ence stars must be taken into account for the most accurate cor¬ 
relation of the FGS data with the radio observations. However, 
the proper motions provided by PPMXL have large errors. As is 
standard in FGS astrometry (e.g. Nelan & Bond 2013, Benedict 
et al. 2007), the PPMXL proper motions are input as “observa¬ 
tions with errors” into the model that combines the FGS astrom¬ 
etry from the eight epochs on to the master plate. The model, 
which applies a 4-parameter plate solution, outputs the best fit¬ 
ting proper motion for each reference star, based upon the FGS 
measurements but constrained by the PPMXL input values. The 
proper motions for HD 93129Aa and HD 93129Ah were solved 
for directly without any external constraints. Table|^lists the ref- 
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erence star identifier, PPMXL coordinates and proper motions 
with the catalog errors, and the resulting FGS tangent plane co¬ 
ordinates and proper motions and errors. The positional errors 
for RA, Dec are in mas, while the proper motion errors are in 
mas yr The FGS ^,7/ coordinates are in arcsecs, and their er¬ 
rors are in mas. The FGS proper motions are now absolute rather 
than relative. We note that PPMXL does not contain reference 
star R2. For this star we used the UCAC4 (Fourth US Naval Ob¬ 
servatory CCD Astrograph Catalog) coordinates but solved for 
the proper motion using only FGS data. 

Next step is to find the RA,Dec of the zero point of the 
FGS tangential plane. The accurate FGS proper motions are used 
to update the PPMXL coordinates to epoch 2009, which is the 
epoch of the FGS master plate. Using the algorithms outlined 
in [Smart ] ( |1960| l, (1) initial values for the RA,Dec of the FGS 
tangent point are given, for which the Guide Star Catalog v2.3 
coordinates of HD 93129A are used, (2) the PPXML (RA,Dec) 
of each reference star are converted to Cartesian coordinates in 
this tangential plane, which are treated as “observations with er¬ 
rors”, (3) the PPXML Cartesian coordinates are used to set the 
scale and orientation of the FGS frame to (East, North). (4) An 
iterative process is applied to find a global solution that provides 
the best-fit of the observed FGS and PPMXL coordinates in 
the tangential plane. The solution creates a hybrid FGS-PPMXL 
catalog and computes the RA,Dec of the tangent point. Once 
convergence is achieved, it is a straightforward process to com¬ 
pute the (RA,Dec) of any point in the FGS frame. We used the 
GAUSSFIT pefferys et al.|1987| l least-squares program to build 
the models to carry out these calculations. 

We compute the RA,Dec of the HD 93129Aa and Ab com¬ 
ponents, adjusted for proper motion to the Aug 6, 2008 epoch of 
the radio observations, equinox J2000; 

a,d(Aa) = 10''43”57.455^-59°32'51.36", 
a,d(Ab) = 10M3”57.456^-59°32'5L33". 

The uncertainty in each coordinate is +27 mas, dominated by 
the PPMXL catalog errors in the position and proper motion of 
the seven FGS reference stars surrounding HD93129A. How¬ 
ever the position of Aa relative to Ab is accurate to about 1 mas. 
From these positions, the separation and position angle of Aab 
in Aug 2008 was 36+1 mas and PA - 12+ 1°. 

The proper motions derived from the FGS data are jUq., /i^(Aa) 
= -8.4,2.6 masyr^' and jUQ.,jUa(Ab) = -9.0,0.0 masyr''. Note 
that the proper motion of the binary’s components include any 
putative orbital motion along with the system motion. The dif¬ 
ference from the reference star proper motions (Table [^, partic¬ 
ularly in RA, could reflect internal motion within the core of Tr 
14, while the reference stars lie outside the core. 

3.2. Orbit estimation 

The relative motion of the components Aa and Ab hints at 
HD 93129A being a gravitationally bound system. However, the 
relative positions shown in Figure 1 cover a rather small part of 
the orbit and thus are inadequate to perform a standard determi¬ 
nation of orbital elements. No definitive radial velocity variation 
has been reported, which is not surprising given the wide sep¬ 
aration and apparent near linear trajectory of the stars. Never¬ 
theless, a first approximation of an orbital fit was attempted. An 
algorithm that minimizes the weighted square distance between 
the measured data points and the fitted orbit (Casco & Vila, pri¬ 
vate communication) was used to determine a preliminary set 
of orbital parameters. Our results suggest an orbital period of 
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(5RA (mas) 

Fig. 4. A preliminary fit for the orbit (triangles). The primary Aa is at 
the origin, The ditferent positions of the secondary Ah between 1996 
and 2013 are represented with crosses, and are the same as in Fig.[T] 
The figure shows the projection of the preliminary orbit in the plane of 
the sky. 


the order of 200 years, an eccentricity larger than 0.9, and semi¬ 
major axes of about 93 and 37 mas respectively for components 
a and b. Our best-fit solution points also to a periapsis argument 
of about 220 degrees and an inclination of about 103 degrees. An 
attempt to fix the inclination value close to zero yielded a much 
poorer fit (larger weighted root mean square rms). According to 
this preliminary solution, the next periastron passage is expected 
to take place in 2024. It should be emphasized that such a pre¬ 
liminary solution was calculated to give a rough idea of the orbit, 
even though the present astrometry does not allow for a convinc¬ 
ing characterization of the orbit accompanied by relevant error 
bars on the estimated parameters. This calculation has however 
the merit to provide a preliminary basis to organize future obser¬ 
vation campaigns dedicated to HD 93129A. 



Fig. 5. Image of the HD93129Aab system, from 2007 LB A data (in 
contours) and the location of stars Aa (southern) and Ah (northern) 
marked with crosses. The contours start at 0.3 mjy bearn^* (3cr) and in¬ 
crease by factors of2F2.The synthesized beam (0.2"x0.05", FA = 31°) 
is shown at the bottom right corner. North is up and East is to the left. 



Fig. 6. The uv coverage of project V19 IB in Aug 2008. 


4. Radio observations 

4.1. Observations of June 22, 2007 

We observed HD93129A on June 22, 2007 (MJD 54273) be¬ 
tween 03:00 and 13:00 UTC at 2.3 GHz as part of eVLBI ex¬ 
periment vtllD3 that used Parkes, Mopra, and Australia Tele¬ 
scope Compact Array (ATCA) to search for suitable phase cal¬ 
ibrators near the target (|Benaglia et al.|2010|l. The sources PKS 
0637-752 (J0635-7516) and J1047-6217 were used as flux 
and phase calibrators, respectively. The total time observing 
HD93129A was 3 hours. 

Standard VLBI calibration procedures were used (see also 
Sect. |4.2| l and the source imaged with a synthesized beam of 
0.2"x0.05", PA = 30.9°. The source has a measured flux density 
of ~ 3 mJy. See more details in|Benaglia et al.|([2010|l and Fig.|^ 


4.2. Observations of August 6, 2008 

Further observations were conducted on Aug 6, 2008 
(MJD 54684) from 01:00 to 12:00 UTC with the LBA at an ob¬ 
serving frequency of 2.3 GHz, and using five antennas: ATCA 
(in tied-array mode, i.e., working as a single station), Parkes, 
Mopra, Ceduna, and Hobart (experiment V191B). The baseline 
lengths range from 100 to 1700 km, providing an angular resolu¬ 
tion of ~ 15 mas at that frequency. The data were recorded at 256 
Mbps provided by 4 sub-bands of 16-MHz bandwidth at each of 
the two circular polarizations. The data were correlated with the 
Distributed FX (DiFX) software correlator located at Curtin Uni¬ 
versity of Technology (WA, Australia), using 128 channels per 
sub-band, each band 125 kHz wide, and using 2-second integra¬ 
tions. The data from the ATCA antennas were also correlated 
as an independent interferometer to obtain the total flux den¬ 
sity from a low-resolution observation. The time on HD 93129A 
amounted to 3.2 hours. The uv coverage is shown in Fig.|^ 
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Table 3. Radio flux density of HD 93129A at different frequencies ob¬ 
tained with ATCA and LBA data. We include the project code and the 
flux density calibrator name with its flux value. 


Array/project 
and Date 

y 

[GHz] 

Sv 

[mJy] 

Flux cal 

name 

^ y,cal 

[Jy] 

ATCA/C678 

2003-01-28 

4.8 

4.1 +0.4" 

1934-638 

2.84 

2003-01-28 

8.6 

2.0 + 0.2" 

1934-638 

5.83 

2003-12-20 

1.4 

9.4 + 0.9" 

1934-638 

14.98 

2003-12-20 

2.4 

7.8 + 0.4" 

1934-638 

11.59 

2004-05-05 

17.8 

1.8 + 0.15" 

Mars 


2004-05-05 

24.5 

1.5 + 0.35" 

Mars 


ATCA/V191B 

2008-08-06 

2.3 

7.5+0.11 

0637-752 

5.32 

ATCA/C1726 

2009-01-18 

4.8 

5.6+0.3 

1934-638 

5.83 

2009-01-18 

8.6 

2.9+0.3 

1934-638 

2.86 

LBA/V191B 

2008-08-06 

2.3 

2.9+0.51*’ 

0637-752 



References, a: Benaglia et al. 2006; b: Derived from image portrayed 
in Fig. [fright. 


The data reduction was performed using the Astronomi¬ 
cal Image Processing System (AlPSjj^ software package. Bad 
data were flagged using telescope a-priori flag files and infor¬ 
mation provided in the observing logs. The amplitude calibra¬ 
tion was performed using the antenna system temperatures and 
phase solutions on the calibrators were obtained using the AIPS 
task FRiNG. The compact source J1047-6217 was used as the 
phase calibrator. The fringe finders 0637-752 (J0635-7516), 
1549-790, J1051-6518, J1023-6646 and the phase calibrator 
J1047-6217 were used for the bandpass calibration. First, an 
accurate model of the calibrator was obtained after several self¬ 
calibration and imaging cycles. The calibrator, slightly resolved 
on the longest baselines, had a flux density of 0.99 + 0.02 Jy. 
Finally, the phase solutions and the amplitude scale were trans¬ 
ferred to the target source, and images produced by deconvolving 
the visibility data. 

The data from ATCA were correlated in standard interferom¬ 
eter mode. The resulting fluxes for the phase calibrator J1047- 
6217 and HD 93129A were 1.300+0.002 Jy and 7.5+0.1 mJy re¬ 
spectively. The flux and uncertainty for the phase calibrator was 
derived by fitting a point source to the image, using all baselines. 


5. Results 

5. 1. Radio flux density 

In Fig. we show the radio image from LBA 2008 observations 
of HD 93129A at 2.3 GHz. Two final images are presented; one 
at maximum angular resolution (Fig. |7]-left: high-resolution im¬ 
age), and a second one using tapering that increases the sensitiv¬ 
ity to larger scale emission (Fig. [7]-right: low-resolution image). 
The high-resolution image has an rms of 0.21 mJy beam ', and 
shows a “comma-shaped” extended emission with a flux density 
of 1.5 + 0.5 mJy. The value quoted for the error is the difference 
in the integrated fluxes above 2 and 3cr, where cr is the measured 
rms of the image. In the low-resolution image the rms noise is 
lower and more flux is recovered from the source, as expected. 
The rms of the image is 0.13 mJy beam \ and a gaussian fit to 
the source gives an integrated flux density detected by the LBA 
of 2.9 + 0.5 mJy, with a peak intensity of 1.8 + 0.2 mJy beam^’. 

We note that the amplitude calibration of the data, which is 
based on the system temperature of the antennas is limited by 
the uncertainty on the antennas’ response and gain curves. The 
flux density measured for the phase calibrator is ~ 30% lower 
than its unresolved flux density measured with ATCA between 
Feb 2007 and Nov 2008; 1.46 + 0.01 Jy to 1.42 + 0.02. This 
indicates that calibrator is partially resolved by the LBA or it is 
variable, or both. We also note that the LBA-measured flux of 
HD 93129A is ~35% of the total flux measured in simultaneous 
Aug 2008 ATCA data at 2.3 GHz. We conclude that the rest of 
the flux from HD 93129A is resolved out by the LBA. 

5.2. Radio morphology and astrometry 

From the 2008 LBA image the angle subtended by the emitting 
source as seen from the secondary star (so-called opening angle 
2flw) can be estimated, leading to half-opening angle 0* ~ 65°. 
This is related to the wind-momentum rate ratio rjhy 0^ = 1 80° x 
7/(1 + rj), thus implying rj ~ 0.6. 

The radio emission at maximum angular resolution presents 
two maxima. The position of the centroid of the two fitted gaus- 
sians is ac = 10;43;57.456, 6c = -59;32;5L34, J2000. 

A systematic uncertainty of 1.1 mas in Right Ascension and 
2.8 mas in Declination comes from the position error of the 
phase reference calibrator in the ICRF2 (International Celestial 
Reference Frame, second realization). Therefore, the absolute 
uncertainty of the position of the radio emission is about 3 mas in 
the ICRF. In Fig. the crosses mark the positions of HD 93129 
Aa and Ah, with the size of the crosses indicating the absolute 
uncertainty, although the relative separation of the components 
has an uncertainty of ~ 1 mas (see Sect.|^. 


4.3. ATCA observations 


6. Discussion 


The Australia Telescope Online Archive (ATOA) was searched 
for observations of HD 93129A. In addition to project C678 (Be¬ 
naglia et al. 2006) and the observations associated with the LBA 
projects described in the previous Section, project C1726 tar¬ 
geted HD 93129A for 3.3 h at 4.8 and 8.6 GHz. Tablej^summa- 
rizes the dates and observing frequencies for these observations. 
In all observations with the ATCA, HD 93129A is not resolved 
by ATCA and appears as a point source. 


6.1. Wind-momentum rates 


The radio source appears in the LBA observations as a bow¬ 
shaped region, slightly curved around the component Ah, and 
reminiscent of the WCRs resolved in WR 140 ([Dougherty et al. 


2005|l, WR 146 ( [O’Connor etal.|2005| l, Cyg OB2 #5 ( [Ortiz-Leon! 


et al. 201 l|l, as examples. If we assume that the radio source is 


centered on the stagnation point of the WCR, the wind momen¬ 
tum rates ratio rj can be expressed as; 




^ http://www.aips.nrao.edu 
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Fig. 7. The LBA observations of HD93129A at 2.3 GHz. a) Left: At highest angular resolution with an image rms of 0.21 mjy beam^' and 
contour levels shown for -0.6, 0.6 (3 it), 0.7, 0.8 and 0.9 mJy bearn^*. b) Right: a tapered image that increases the weight of the shorter baselines 
and increases sensitivity to larger scale emission; the image rms is 0.13 mJy bearn^* and the contour levels are -0.4, 0.4 (3 (t), 0.6, 1, 1.4 and 1.8 
mJy beam^'. The synthesized beams are 15mas x llmas and 31mas x 23mas respectively, shown in the bottom right-hand comer of the images. 
The crosses mark the positions of the system components Aa (south) and Ah (north), at the epoch of the radio observations based on the position 
derived in Sec.[^ North is up and East is to the left. 


where R, is the distance between the WCR and star i ( |Usov| 
1992| l. Note that the value of j] is independent of orbit inclina¬ 
tion. 


From the relative position of Aa, Ah and the centroid of the 
radio emission in Figure 7, we estimate RbIRa ~ 0.7. The stellar 
winds of the two stars are assumed to have reached their termi¬ 
nal velocities before colliding, considering the wide orbit, and 
we further assume that the terminal velocities of both compo¬ 
nents are similar (see the values quoted for instance in Table 1 
of Muijres et al. 2012), then 77 reduces to a mass-loss rate ratio, 
leading to Mi, 0.5Ma¬ 


in this work, we hesitate to use the shape of the WCR to es¬ 
timate 77 . As Pittard & Dougherty ( 2006) l (section 4.3.3) show 
for WR 140, the emission in the WCR detected in VLBI is close 
to the stagnation point and the opening angle has not reached 
the asymptotic value. They also noted the challenge of identi¬ 
fying the stagnation point relative to the stellar components due 
to opacity effects. In the best case, we can perhaps assume the 
shape of the contact discontinuity (CD) between the two stellar 
winds (jCanto et al.|1996| l is represented by the radio emission. In 
Figure)^ the positions of the CD corresponding to different val¬ 
ues of 77 are plotted relative to the emission detected in the LBA 
observation. The distance and angle between stars are known, so 
the only free parameters are the offsets with respect to the radio 
image and 77 . The 77 values that best match the radio source are 
between ~ 0.4 and ~ 0.6. 


6.2. On the binary nature of HD 93129A 


Do HD 93129Aa and Ah form a gravitationally bound binary 
system? The system is being monitored in the OWN Survey 
project ( Barba et al.||2010 l. A radial velocity (RV) curve from 
preliminary data shows slight variations though no period has 
yet been determined. The putative low amplitude RV variations 
may be due to a number of factors, including insufficient data 
over a potential orbital period, a very low mass companion or 



Fig. 8. Image of the HD 93129Aab system at 2.3 GHz, showing different 
77 values. The black dots mark the position of the stars. The synthesized 
beam is shown in the bottom right comer. 


unfavorable orbital inclination (Barba & Gamen, private com¬ 
munication). 

An argument in support of HD 93129Aab being a bound sys¬ 
tem is the small angular separation and the change in the rate at 
which the two stars are approaching one another. In the simplest 
scenario, if they are in the plane of the sky, with a separation of 
26.5 mas (epoch 2013.09 from Sana et al. 2014) and at a distance 
of 2.5 kpc, they have a physical separation of only 66 AU, which 
implies a bound system for such massive stars. However, this is 
strongly dependent on the inclination being in the plane of the 
sky. More compelling is the increase in the rate at which the two 
stars are approaching one another. Based upon FGS and ACS 
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data, between 1996 and 2009 the separation had been closing 
at approximately 2.4 mas yr *, but VLTI data shows this rate to 
have accelerated to 4.2 mas yr"' between 2011.19 and 2013.09 
( |Sana et al.|2014| l, consistent with a binary system approaching 
periastron. The relative positions in Fig. [T] suggest either an or¬ 
bital inclination close to 90 degrees, or a nearly face-on orbit of 
high eccentricity. The results of the very preliminary orbit fits 
given in Sect.[^favor the former hypothesis. 


6.3. On the mass-loss rates of the HD93129A system 

HD 93129A has been the tar^ of different studies a imed at de¬ 
riving the mass-loss rate (MF 


Taresch et al. 


(19971 estimated 


M of HD 93129A assuming it was a single star, by htting the 
Ha profile and also using a more complete spectral synthesis 
approach in the ultraviolet. These two approaches led to con¬ 
sistent values of 1.8 x 10“^ and 2.1 x 10 “^ Mq yr“\ respectively. 
The authors considered that these values should be reduced by 
a factor 2^^^, a reduction of about 40%, if HD93129A con¬ 
sisted of two very similar components. Using Ha line diag¬ 
nostics exclusively, Repolust et al. (2004) derived a value of 
2.4 X 10“^ Mg yr“' for the mass-loss rate, noting that the results 
could strongly suffer from contamination by a companion star. 

[Benaglia & Koribalski| ( |2004| l presented flux density mea¬ 
surements of the HD93129A system from ATCA data at two 
frequencies, 4.8 and 8.64 GHz, and derived a spectral index 
a = -1.2 + 0.3. Their interpretation was that both thermal and 
non-thermal emission were contributing to the radio emission. 
[Wright & Barlow] ( |1975| ) showed that the mass-loss rate of a 
star with a thermal wind can be expressed in terms of the radio 
flux density. Using the Wright & Barlow (19751 formulae and 
if /t is the fraction of thermal emission at a given frequency, 
it is possible to express the mass-loss rate in terms of fj as 
M = X 7.2 X 10-5 Mg yr-' at v = 8.64 GHz for a stellar 

distance of 2.5 kpc. This leads to a mass-loss rate of ^ 5 x 10 ^ 
Mgyr-'if/T = 0.5. 

At the epoch of the observations presented by |Benaglia~& 
|Koribalski| ( |2004| l, the stars Aa and Ab were ~ 50 mas apart in 
the plane of the sky (see Sect. 2). The angular resolution of the 
ATCA data at ~ 1 arcsec precludes separation of the contribu¬ 
tions to the total flux of both stellar winds and a colliding-wind 
region. An estimate of the thermal emission from the two stellar 


winds can be obtained from the radio observations in Benaglia et 
jak] ( |2006) from ATCA at 17.8 and 24.5 GHz. At these frequen- 
cies, non-thermal emission is negligible and can be disregarded 
to zeroth order. The total flux is then the sum of the thermal emis¬ 
sion from the two stellar winds, and characterized with a spectral 
index of h- 0.6 ( Wright & Barlow|1975 1. Consequently, the ther¬ 
mal contribution to the flux at lower frequencies can be derived 
by extrapolation and assuming no variations between 2003 (8.64 
GHz data) and 2004 (17.8 GHz data). This leads to a thermal 
contribution at 8.64 GHz of 1.2 mJy, compared with a total flux 
of 2.00+0.15 mJy (see Table[^, an /x ~ 0.6, and M - 4.9 x 10^5 
Mgyr-'. 

The main conclusion here is that mass-loss rate determina¬ 
tions are challenging when dealing with a binary system. Con¬ 
sidering the high fraction of binaries among O-type stars (see 
e.g. jSana et al.|2014] and references therein), many observational 
determinations of M should certainly be viewed with caution, in¬ 
cluding those of HD 93129A. 


® For an overview of these mass-loss rate determination techniques, in¬ 
cluding a comprehensive discussion of their advantages and limitations. 


we refer to the review by Puls et al. j2008|. 
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Table 4. Mass loss rates values derived for HD 93129A. 


Reference 

Method 

Value xlO^ 
(Mg yr-') 

Tar97 

Ha profile 

1.8 

Tar97 

Ultraviolet lines 

2.08 

Rep04 

Ha profile 

2.36 

BK04 

8.64GHz continuum flux 

(/t) 3/4 7.2t 

This work 

Separating T and NT fluxes 

4.9tt 


References. Tar97: [Taresch et al.|(|1997|(; Rep04:|Repolust et al.[(|2004|l; 
BK04: [Benaglia & Koribalski[^2004| >; j: (fy) e [0,1]; ft: see text. 


6.4. Flux density variability of HD 93129A 


The flux density measurements of HD 93129A derived from 
ATCA data (Table cover frequencies from 1.4 to 24.5 GHz. 
Observations at 2.4,4.8 and 8.6 GHz were obtained twice, at dif¬ 
ferent epochs. The two measurements at 2.4 GHz were observed 
in 2003 and 2008, and the fluxes agree within the uncertainties. 

The fluxes at 4.8 GHz (from 2003 and 2009) differ by ~ 
37%, increasing with time. Similarly, the 8.6 GHz fluxes at the 
same epochs differ by ~ 45%. In all the archival observations, 
relative errors are 10% or less. We checked for calibration errors 
in the individual datasets. For all of them, the flux calibrator was 
1934-638, from which the flux scale of phase calibrators and 
source were bootstrapped. 

The data indicate that the flux increased by a significant frac¬ 
tion compared to the uncertainties between 2003 and 2009. Such 
a trend is indeed anticipated if the system is approaching perias¬ 
tron passage c.f. WR 140 ( [Dougherty et al.|2005[ ). 

A synchrotron component that dominates over the thermal 
emission from the stellar winds of the two stars can explain the 
measured flux densities. We have modeled the data with a power- 
law spectrum with a low-frequency exponential cutoff. We have 
computed two possible models. The first model fits the data from 
January 2003 - December 2003, and the 2004 data at the higher 
frequencies (Fig. 9; solid dots). It is assumed that the flux does 
not vary significantly within this range of time, given the sep¬ 
aration of Aa to Ab in that period. These data are well-fit by a 
power-law, with a turnover around 1.4 GHz. However, the data 
from 2008 and 2009 are modeled by a similar power-law spec¬ 
trum but with a turnover at a slightly higher frequency. An addi¬ 
tional thermal component is required to fit the data at 20 GHz. 
The spectral indices from the fits to the 2003-4 and 2008-9 data 
are, respectively, -1.03 + 0.09 and -1.21 + 0.03. 


Up to now, the results from observations at high-energy spec¬ 
tral bands remain inconclusive in detecting emission from a 


WCR. In soft X-rays (below 10 keV), the work by Cohen et al. 
P011[ l concluded that the X-ray spectrum is dominated by emis¬ 
sion from the individual stellar winds in the system. On the other 
hand, Gagne et al. (20111 reported a ratio between the X-ray and 
bolometric luminosities (Lx/Cboi) of 1.3 x 10 “^ xhis value sug¬ 
gests that only a small contribution of the X-rays come from the 
colliding-wind region. It is not unexpected to see colliding-wind 
binaries made of very early-type stars with no bright X-ray emis¬ 
sion in excess of the canonical Lx/Cboi value of ~ 10^’. For very 
early-type stars with strong and dense winds such a the members 
of HD93129A, the individual X-ray emission should be some¬ 
what fainter than suggested by the canonical ratio (Owocki et 
al. 2013; De Becker 2013). A low emission measure of a WCR 
due to a large stellar separation on top of such an undeiTuminous 
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Fig. 9. Spectra of HD 93129A obtained from ATCA data at different 
epochs (see labels). The data from 2003 have been fitted with two dif¬ 
ferent models: a Free-Free Absorption (FFA, dotted line), and a FFA 
plus a thermal emission also considering the 2004 data (solid line). The 
2008-2009 data have been fitted with a FFA model (dashed line); see 
also Pittard & Dougherty (2006). 


wind emission leads therefore to a not so high overall Lx/^-boi 
ratio. 

HD93129A is rather close to the Fermi source IFGL 
J1045.2-5942 (in the 0.1-100 GeV energy range). This source 
is associated with tj Car, the only PACWB identified on the basis 
of y-ray emission. It seems very unlikely that HD 93129A con¬ 


tributes to the emission detected by Fermi, as clarified by Abdo 
|et al.| ( |2010[ ). At even higher energies, no TeV source has been 
found at the position of HD93129A in the TevCat catalogue 
(http;//tevcat.uchicago edu/), even when considering newly an¬ 
nounced sources and source candidates. We conclude that there 
is no known y-ray source associated with HD93129A. Obser¬ 
vations with the forthcoming Cherenkov Telescope Array, with 
unprecedented sensitivity and source location techniques, will 
surely provide useful constraints to better understanding of the 
high energy phenomena in these objects. 


7. Summary and conclusions 

An LBA observation of the massive binary HD 93129A detected 
an extended and arguably curved radio emission component with 
a flux density of 2.9 mJy at 2.3 GHz between the two compo¬ 
nents Aa and Ah. Following a detailed analysis of recent high 
angular resolution astrometry of the system, we provide com¬ 
pelling evidence that the radio emission component is coincident 
with the expected position of a wind interaction region between 
components Aa and Ah, and suggest a wind-momentum ratio of 
~ 0.5. 

Historic ATCA observations of the HD 93129A system, from 
1 to 25 GHz, re-reduced in a uniform way, show that the flux in¬ 
creased between 2003-4 and 2008-9, both epochs being well fit 
by a power-law spectrum with a steep spectral index. Similar 
increases in flux have been seen in other systems as they ap¬ 
proach periastron. Thus, the results presented in this paper lend 
significant additional support to the idea that wind-collision re¬ 
gions are the sites where relativistic particles are accelerated in 
particle-accelerating colliding-wind binaries. 

Very Long Baseline Interferormeter (VLBI) observations 
of particle-accelerating colliding-wind binaries help to quan¬ 
tify specifically the properties of the non-thermal radio emis¬ 
sion. Such measurements are crucial for models aimed at re¬ 
producing the particle acceleration and non-thermal physics at 


work in these objects. VLBI observations of the HD 93129A 
system across a wide frequency range will allow more accurate 
properties of the spectrum to be derived, and hence the proper¬ 
ties of the relativistic electron populations involved in the syn¬ 
chrotron emission process, and potentially reveal associations 
between flux variations and changes in the properties of the 
wind-collision region. 

In parallel, new ATCA observations, preferentially simulta¬ 
neous across all the observing frequencies are necessary to ob¬ 
tain accurate total fluxes. Repeated observations in the leading 
up to the putative periastron passage will reveal if the flux in¬ 
creases as observed in other PACWB systems at this phase of 
the orbit. 

Acknowledgements. The authors thank an anonymous referee for detailed com¬ 
ments, C. Phillips and A. Tzioumis for their work in a preliminary stage of the 
study, A. Lopez-Sanchez for providing some of the ATCA observations, Rodolfo 
Barba, Roberto Gamen, Ed Fomalont, Jamie Stevens, and very especially G. 
Vila and N. Casco for help with the orbits. P.B. is supported by the ANPCyT 
PICT- 2012/00878 and UNLP G11/115 project. B.M. acknowledge support by 
the Spanish Ministerio de Economia y Competitividad (MINECO) under grants 
AYA2013-47447-C3-i-P and BES-2011-049886. This study was partially sup¬ 
ported by NASA through grant GO-10898 from the Space Telescope Science In¬ 
stitute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. 
The research has made use of the Washington Double Star Catalog maintained 
at the U.S. Naval Observatory. 


References 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, ApJ, 723, 649 
Barba, R. H., Gamen, R., Arias, J. I., et al. 2010, Rev. Mex. Astron. Astrof. Conf. 
Ser. 38, 30 

Benaglia, R & Romero, G.E. 2003, A&A, 399, 112 
Benaglia, P. & Koribalski, B. 2004, A&A, 416, 171 

Benaglia, R, Koribalski, B., & Albacete Colombo, J. F. 2006, RASA, 23, 50 
Benaglia, R, Dougherty, S. M., Phillips, C., Koribalski, B., & Tzioumis, T. 2010, 
Rev. Mex. Astron. Astrofis. Conf. Ser., 38, 41 
Benedict, G. F, McArthur, B. E., Eeast, M. W, et al. 2007, AJ, 133,1810 
Canto, J., Raga, A. C., Wilkin, F. P. 1996, ApJ, 469, 729 
Cohen, D. H., Gagne, M., Leutenegger, M. A., et al. 2011, MNRAS, 415, 3354 
Contreras, M. E., Rodriguez, L. F, Tapia, M., et al. 1997, ApJ, 488, L153 
De Becker, M. 2007, A&ARv, 14, 171 
De Becker, M. 2013, NewA, 25, 7 
De Becker, M. & Raucq, F. 2013, A&A, 558, A28 

Dougherty, S. M., Beasley, A. J., Claussen, M. J., Zauderer, B. A., & Boling- 
broke, N. J. 2005, ApJ, 623, 447 

Dougherty, S. M., & Pittard, J. M. 2006, in Proceedings of the 8th Euro¬ 
pean VLBI Network Symposium, held in Torun (Poland), eds. W. Baan, R. 
Bachiller, R. Booth, et al. 

Eichler, D., & Usov, V. 1993, ApJ, 402, 271 
Gagne, M., Fehon, G., Savoy, M. R., et al. 2011, ApJS, 194, 5 
Hog, E., Eabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 
Jefferys, W. H., Eitzpatrick, M. J., McArthur, B. E. 1987, CM, 41, 39 
Maiz Apellaniz, J., Walbom, N. R., Morrell, et al. 2008, Rev. Mex. Astron. As¬ 
trofis. Conf. Ser., 33, 55 

Mason, B. D., Gies, D. R., Hartkopf, W. I., et al. 1998, AJ, 115, 821 
Meunier, C., Bauer, E. E., Dwarkadas, V. V., et al. 2013, MNRAS, 431, 2453 
Muijres, L. E., Vink, J. S., de Koter, A., Muller, P. E., & Danger, N. 2012, A&A, 
537, A37 

Nelan, E. R, Walborn, N. R., Wallace, D. J., et al. 2004, AJ, 128, 323 
Nelan, E. R, Walborn, N. R., Wallace, D. J., et al. 2010, AJ, 139, 2714 
Nelan, E.P. & Bond, H.E., 2013, ApJ, 773,L26 

Nelan, E.P., 2014, Fine Guidance Sensor Instrument Handbook, Version 22.0 
(Baltimore, MD:STScI) 

O’Connor, E., Dougherty, S., Pittard, J., & Williams, P. 2005, in Proceedings of 
"Massive Stars and High-Energy Emission in OB Associations", a workshop 
of the JENAM 2005, "Distant Worlds", held in Liege (Belgium), July 4-7, 
2005, Eds. G. Rauw, Y. Naze R. Blomme, and E. Gosset, p. 81 
Ortiz-Leon, G. N., Loinard, L., Rodriguez, L. F, Mioduszewski, A. J., & Dzib, 
S. A. 2011, ApJ, 737, 30 

Owocki, S. R, Sundqvist, J. O., Cohen, D. H., & Gayley, K. G. 2013, MNRAS, 
429, 337 

Payne, C. H, 1927, Harvard Bull. 846, 10 

Pittard, J. M., & Dougherty, S. M. 2006, MNRAS, 372, 80 

Puls, J., Vink, J. S. & Najarro, F. 2008, A&ARv, 16, 209 


Article number, page 9 of 


10 














A&A proofs: manuscript no. aa25595-astroph 


Repolust, T., Puls, J., & Herrero, A. 2004, A&A, 415, 349 
Roeser, S., Demleitner, M., Schilbach, E. 2010, AJ, 139, 2440 
Sana, H., Le Bouquin, J.-B., Lacour S., et al., 2014, ApJS, 215, 15 
Smart, W. M., "Text-Book on Spherical Astronomy", 4th Edition 1960, Cam¬ 
bridge University Press 

Sota, A., Mai'z Apellaniz, J., Morell, N. E, et al. 2014, ApJS, 211, 10 
Stevens, I. R., Blondin, J. M. & Pollock, A. M. T., 1992, ApJ, 386, 265 
Taresch, G., Kudritzki, R. R, Hurwitz, M., et al. 1997, A&A, 321, 531 
Usov, V. V. 1992, ApJ, 389, 635 
Walborn, N. R. 1982, AJ, 87, 1300 
Walbom, N. R. 1995, RMxAA, 2, 51 

Walborn, N. R., Howarth, I. D., Lennon, D. J., et al. 2002, AJ, 123, 2754 
Williams, P. M., Dougherty, S. M., Davis, R. J., et al. 1997, MNRAS, 289, 1 
Wright, A. E. & Barlow, M. J. 1975, MNRAS, 170, 41 


Article number, page 10 of 


10 



